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(57) Abstract 

The present invention relates to a method for determining characteristic electrical properties of semi-conducting materials wherein 
the tune/frequency dependent electrical impedance or admittance of the material is measured. The invention also relates to an apparatus 
for carrying out the method. A number of bulk and surface parameters characterize the electrical properties of a given piece of material 
These parameters include the dielectric constant e of the material, the difference Ap<*> in the chemical potential of the bulk of a material 
and the chemical potential of its surface and/or of metal electrode - material surface interface, the density of the majority and minority 
electrical mobile charge carriers N and Nmm, respectively, in the bulk of the material, the electrical mobility n of the majority electrical 
mobile charges in the bulk of the material and the electrical mobility of minority mobile charge carriers, the surface and bulk emission 
and capture rates ESL and CJL, respectively, for mobile positive and negative charge carriers characterizing the effect of surface and bulk 
localized states within the band gap, when they are present, on the electrical transport 



FOR THE PURPOSES OF INFORMATION ONLY 



Codes used to identify States party to the PCT on the front pages of pamphlets publishing international 
applications under the PCT. 



AT 


Austria 


GB 


United Kingdom 


MR 


bCsubauia 


AD 


Australia 


GE 


Georgia 


MW 


Malawi 


BB 


Barbados 


GN 


Guinea 


NE 


Niger 


BE 


Belgium 


GR 


Greece 


St 


Netherlands 


BF 


OUffctUSA FSSO 


eu 


Hungary 


NO 


Norway 


BG 


Bulgaria 


IE . 


Ireland 


NZ 


New Zealand 


BJ 


Boom 


FT 


Italy 


PL 


Poland 


BR 


Brazil 


JP 


Japan 


FT 


Portogal 


BY 


Belarus 


KE 


Kenya 


RO 




CA 


Canada 


KG 


Xytgystan 


RIT 


Rmgrisu Pfdffattoo 


CF 


Central African Repubflc 


KP 


Democratic People's Republic 


SD 


Sudan 


CG 


Congo 




of Kojca 


SE 


Sweden 


CH 


Switzerland 


KR 


Republic of Korea 


SI 


Stowsia 


a 


Cotedlvobe 


KZ 




SK 


Slovakia 


CM 


Cameroon 


U 




SN 


Spnqgai 


CN 


China 


LK 


Sri Lanka 


ID 


Chad 


CS 


Czechoslovakia 


UJ 


Lmccobourg 


TG 


Togo 


CZ 


Czech Republic 


LV 


Utvia 


ti 


Tajikistan 


DE 


Gammy 


BflC 


Monaco 


IT 


THnldad and Tobago 


DK 


Denmark 


MD 


Republic of Moldova 


UA 


Ukraine 


ES 


Spain 


MG 


Madagascar 


US 


UnUed States of America 


FI 




ML 


Mali 


uz 


Uzbekistan 


FR 


Vance 


MN 


Mongo&a 


VN 


Viet Nam 


GA 


Gabon 











WO 94/14078 



PCI7DK93/00410 



1 

METHOD AND APPARATUS FOR DETERMINING CHARACTERIS- 
TIC ELECTRICAL PROPERTIES OF SEMI -CONDUCT I NG 
MATERIALS - 

The present invention relates to a method for determining 
5 characteristic electrical properties of semi -conducting materials 
wherein the time/ frequency dependent electrical impedance or 
admittance of the material is measured. The invention also relates 
to an apparatus for carrying out the method. 

10 A number of bulk and surface parameters characterize the electrical 
properties of a given piece of material. These parameters include 
the dielectric constant t of the material, 

rh 

the difference Am in the chemical potential of the bulk of a 
material and the chemical potential of its surface and/or of metal 

15 electrode - material surface interface, 

the density of the majority and minority electrical mobile charge 
carriers N and N min> respectively, in the bulk of the material, 
the electrical mobility fi of the majority electrical mobile charges 
in the bulk of the material and 

20 the electrical mobility p m j n of minority mobile charge carriers, 

the surface and bulk emission and capture rates E.R and C.R., 
respectively, for mobile positive and negative charge carriers 
characterizing the effect of surface and bulk localized states 
within the band gap, when they are present, on the electrical 

25 transport. 

Depending on the experimental arrangement the present invention can 
be used either to characterize the whole piece of material (integral 
characterization), or alternatively to map the variations of the 

30 above mentioned electrical material parameters, i.e. some of them, 
all of them and/or various combinations of them such as the specific 
resistivity p = 1/N.e.p as a function of the position of the 
measuring head and measuring point on the surface of the material be 
that a piece of semi-conductor, insulator or a semi -conducting 

35 wafer. 

There exists a number of experimental methods today that are being 
used for the experimental determination of the above mentioned 

rh 

electrical material parameters Ln , £, N and fi, in semi-conductors 
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and insulators. Table 1 illustrates to the best knowledge of the 
inventors the range of the experimental methods (the list not being 
exhaustive) that are relevant to the present invention. 

The literature about each of these experimental methods is vast. The 
references (1), (2) and (3) relate to Standard Electrical Impedance 
Experimental method and have been chosen as sufficient background 
reference since the present invention relates directly to this 
particular experimental method. 

The known experimental methods used today for the determination of 
the parameters suffer from a number of disadvantages the main of 
which can be summarized as follows: 



15 a) None of the existing experimental methods is capable of 
determining all of the parameters simultaneously and 
independently in a single experiment. 

b) A combination of two different experimental methods is often 
required in order to be able to determine some of the parameters 

20 in question, for example a combination of the Hall effect and the 

d.c. four probe resistivity measurements is required if both n 
and N are to be obtained. 

c) Determination of some of the parameters in question is often very 

ch 

indirect. The determination of Ajr by Photoemission, for 

25 example, is measured in two steps. First the work function of the 

un-contacted surface of the material is measured and then the 

measurement of the work function of the contacting metal is 

performed. The difference between the two work functions is then 

ch 

assumed to be a measure of the A/r when the material in question 
30 is contacted by the given metal. 

d) The experimental methods in question are often very expensive to 
establish (PEE,FTIR,FTIPL), requiring either an ultra high vacuum 
environment (CP,TEE,PEE), liquid Helium temperatures (FTIR,FTIPL) 
or magnetic fields (HE). 

35 e) Because of the very nature of the experimental arrangements the 
existing methods can often be considered only as qualitative in 
relation to the determination of the numerical values of the 
parameters in question, 
f) All the electrical experimental methods named in Table 1 
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(Standard EIS, 4PDCR, HE and TOF) either neglect totally or treat 
inadequately the effects of the electrical contact regions, the 
effects of the real surfaces of the material under investigation 
(surfaces that are part of the electrical measuring circuit), the 
5 effects of the metal electrode - material surface regions and the 

effects of the depletion regions within the material under 
investigation on the calculated values of the parameters in 
question. This neglect renders them unusable for the 
characterization of higher resistivity materials (/> z 100 Q»cm) 
10 and introduces relatively large errors even for characterization 

of materials with resistivities in the region 1.0 < p < 100 G*cm 
(averaging of the effect of electrical contacts, depletion 
regions and bulk). 

15 Among other experimental methods (not displayed in Table 1) used 
today for the characterization of semiconductors there are two that 
are relevant as far as the present invention is concerned (see 
EP-A-0477795, DE-A-3103611, US-A-5103183 and US-A-4028207) and which 
should therefore be mentioned. 

20 

These are the C-V and DLTS experimental methods and variations 
thereof aiming at the determination of either the doping density 
and/or doping density profiles (C-V) or at the characterization of 
the deep levels in semiconductors (DLTS). 

25 

The measured quantity In both methods is an appropriate capacitance 
and its variation with external d.c. voltage bias (C-V) or 
temperature (DLTS). 

30 The major drawback of both is that the measuring frequency at which 
the wanted capacitance is measured is chosen arbitrarily. This fact 
often renders them unusable. 

The proper analysis of the electrical responses embodied in the 
35 present invention is essential in deciding on their validity. This 
prior art in comparison with the present invention is discussed in 
more detail at the end of this specification. 

The above drawbacks of the prior art are eliminated by the method 
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according to the present invention which is characterized by the 
steps of establishing a measurable difference in chemical potential 
between the inner and surface of the material or between the inner 
of the material and the metal electrode - material surface interface 
5 layer, and of determining from the measured electrical impedance or 
admittance one or more of the following electrical parameters of the 
material: 

ch 

the difference A/r in chemical potential, 
the dielectric constant e, 
10 the density N of majority mobile charge carriers, 

the density N m ^ n of minority mobile charge carriers, 
the electrical mobility p of majority mobile charge carriers, 
the electrical mobility /i^ n of minority mobile charge carriers, 
the emission and capture rates E.R. and C.R., respectively, for 
15 mobile positive and negative charge carriers characterizing the 
effect of surface and bulk localized states within the band gap Eg, 
when they are present, on the electrical transport and the measured 
electrical Impedance. 

One of the essential aspects of the present invention is the reali- 
zation that the electrical characteristics of the surfaces of the 
material and/or of metal electrode-surface interfaces, and the 
depletion regions within the material and of the bulk region of the 

material are all interrelated in a calculable fashion that depends 

eh 

on the values of the electrical material parameters Ajr , c, N, p, 
E.R. and C.R. that characterize the material electrically. Therefore 
by including all these regions in the analysis of the measurement of 
the electrical response of the material using for example Standard 
EIS experimental method that leads to the determination of the 
complex electrical impedance Z (w) of the material, the electrical 
material parameters of the material can be determined unambiguously. 

In other words the electrical characteristics of the electrical 
contacts, the depletion regions and the bulk of the material are all 
35 determined by a few fundamental physical parameters that 
characterize the electrical properties of a given piece of material 
- be that semiconductor single crystal ingot, a piece of 
polycrystalline semiconductor, semiconductor wafer or a piece of 
insulating material. 



20 



25 



30 
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Contrary to all the relevant existing electrical experimental 
methods used today for similar type of measurements - methods that 
either neglect, average and/or treat inadequately the effects of 
electrical contacts and depletion regions on the resulting measured 
5 electrical characteristics and which in their experimental arrange- 
ment all aim to achieve "good" ohmic contacts (the effect of the 
electrical contacts and depletion regions being small and therefore 
often neglect able) - the present invention is based precisely on the 
opposite requirement, namely the requirement of establishing non- 
10 ohmic ("bad") electrical contacts, where the electrical 
characteristics of the electrical contacts and depletion regions are 
apparent and are therefore easily measurable. 

The experimentally observable quantity in question that contains all 
15 the necessary information from which the calculation of the 
electrical material parameters is possible is the above mentioned 
complex electrical impedance Z (w). 

To fulfil the requirement of non -ohmic electrical contacts and the 
20 requirement of the possibility of a quantitative determination of 
the electrical material parameters, a number of necessary steps has 
to be taken including the definition of the sample (the piece of the 
material under investigation subjected to a measurement geometry) 
and establishment of the necessary conditions for the existence of 
25 the depletion regions within the sample. 

The present invention therefore relates to a method consisting of a 
number of stages. These are: 

30 1. Sample preparation Definition of the A/i ch for the 



35 



4. Method of analysis 



3. Method of measurement 



2. Experimental arrangement 



sample 

Definition of the sample geometry 
In the given electrical measurement 
circuit 

Definition of the method to measure 
the electrical response of the 
sample 

Calculation of the electrical 
material parameters that 
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characterize the electrical 
properties of the sample. 

The method according to the present invention will be somewhat 
5 different in the two basic applications envisaged. These 
applications are named here as "research" and "on - line quality 
control" configurations, each being characterized by somewhat 
different content of the steps of the present invention. The 
principle of the present invention 1s though the same in both these 
10 configurations. 

The "research" application of the present invention involves more 
detailed studies of the electrical material parameters for example 
as functions of temperature, pressure, heat treatment, electrical 
15 d.c. voltage bias, optical excitation etc. 

This type of studies require a proper treatment of the active 
surface of the material, deposition of the appropriate metal 
electrodes on such prepared surfaces and the experimental 

20 environment (arrangement) that fulfils the above mentioned 
requirements concerning the external parameters, such as 
temperature, pressure etc. If an x,y mapping Is required, this will 
be done destructively by cutting a given piece of material 
(semiconductor wafer) into a number of elements (x,y) and measuring 

25 these sequentially. 

The "on-line quality control" application of the present invention 
on the other hand is Intended for the quality control (integral or 
x,y mapping) of semiconductor wafers. The electrical contacts at a 

30 given x,y position have to be made quickly (non-destructively) and 
after the measurement the measuring head (point) has to move 
reasonably quickly to a new x,y position, so that x,y mapping of the 
wafer of one or more of the electrical material parameters can be 
made reliably and fast and with sufficient throughput (number of 

35 wafers processed in unit of time). This application of the present 
invention therefore involves no or very little sample preparation 
(production wafers) and fast, automated sample exchange mechanism 
together with x,y scanning of the measuring head across the surface 
of the wafer. 
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The apparatus for carrying out the method comprises an external 
electrical circuit for measuring the time or frequency dependent 
electrical impedance or admittance of the material and means for 
electrically contacting the surface or surfaces of the material and 
5 the apparatus is characterized by means for providing a measurable 
difference in chemical potential between the inner of the material 
and its surface or between the inner of the material and the metal 
electrode - material surface interface layer, and means for 
determining from the measured electrical impedance or admittance one 

10 or more of the following electrical parameters of the material: 
the difference Aji ch 1n chemical potential, 
the dielectric constant £, 
the density N of majority mobile charge carriers, 
the density N min of minority mobile charge carriers, 

15 the electrical mobility p of majority mobile charge carriers, 

the electrical mobility /t fflin of minority mobile charge carriers, 
the emission and capture rates E.R and C.R., respectively, of mobile 
positive and negative charge carriers characterizing the effect of 
surface and bulk localized states within the band gap Eg, when they 

20 are present, on the electrical transport and the measured electrical 
impedance. 

The invention will now be described in more details with reference 
to the drawings in which 

25 

Fig. 1 is a schematic diagram of the one-electron energy 

levels in the bulk of the material, at the surface 
of the material and in the metal electrode that are 
relevant to the present Invention. 
30 (Definition of the A/i Ch ) 



Fig. 2a, 2b, 2c are schematic diagrams of possible capacitive 
couplings between the material and the external 
electrical field (external electrical 
35 metal electrodes) that may be utilized to implement 

the present invention. 
2a) "Research" configuration 

("sandwich" capacitive coupling) 
2b) "on-line quality control" configuration 
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("sandwich" capacitive coupling) 
2c) "Research" and/or "on-line quality control" 
configurations 

("planar 0 capacitive coupling) 
5 (Definition of the sample geometry) 



Fig. 2d,2e are schematic diagrams showing the geometrical 

dimensions of the material under investigation 
2d) "Research" configuration 
10 2e) "on-line quality control" configuration 

(Definition of the sample geometry) 



Fig. 3 is a schematic diagram of the experimental 

measuring arrangement which may be utilized to 
15 measure the electrical complex impedance Z (a) of 

the series combination - sample bulk, sample 
depletion regions and sample electrical constant 
regions (sample surface-metal electrode interface). 
(Definition of the measurement method) 

20 

Fig. 4a,4b,4c,4d,4e is a schematic diagram of the passive R,C 
electrical network model that approximates the 
electrical response of the material (bulk, 
depletion and electrical contact regions) to a 
25 small external electrical field input signal. 

(Definition of the approximative R,C model for the 
electrical response function Z s (o>) of the sample) 

Fig. 5 is a schematic graphical representation of the 

30 electrical response function Z (o>), (Y(w)) of the 

material sample modelled by fig. 4d, comprising the 
bulk, the depletion and electrical contact regions 
of the sample. 



35 Fig. 6a, 6b show Example 1; the experimentally observed 

electrical response of the sample 1 (ultra pure 
silicon monocrystal as a function of temperature T. 
(w=2**f) 
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6a) The logarithm of the real part of the electrical 

impedance Z (w) versus logarithm of the measuring 

frequency w, 
6b) The logarithm of C (imaginary part of the 
5 electrical admittance Y(w) divided by (<a) versus 

logarithm of measuring frequency w. 

(Results - example 1) 



Fig. 7a, 7b show Example 2; the experimentally observed 

10 electrical response of the sample 2 (pure, 

polycrystalline silicon, unknown type, unknown N) 
(a=2ff*f) 

7a) The logarithm of the real part of the electrical 
15 impedance Z (*>) versus logarithm of the measuring 

frequency w. 

7b) The logarithm of C(w) (imaginary part of the 
electrical admittance Y(w) divided by u) 
versus logarithm of the measuring frequency a. 
20 (Results - example 2) 

Fig. 8a, 8b, 8c The theoretical analysis of the electrical response 
of the pure, polycrystalline silicon sample 
(Example 2) 

25 8a) A topological model of polycrystalline silicon 

sample consisting of Nj x N 2 grains, each grain 
being characterized by its electrical impedance 
Zj«(u). The shown passive, static R,C electrical 
network representing the total electrical response 

30 includes also metal -semiconductor interface region 

represented by ^h+S' ^m+s elements. 
8b) A simplified physical model of depleted grain "j". 
8c) Static R,C electrical network representing the 
electrical response of a single grain. 



35 



The present invention is not limited to any particular apparatus, 
sample preparation and measurement method but may be carried out by 
a variety of different steps (method stages). 
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It is therefore desired that the invention not be limited to the 
preferred embodiment (figures 1 to 8) and it is intended to cover in 
the appended claims all such variations or modifications as fall 
within the true spirit and scope of the present invention. 

5 

The true spirit and scope of the present invention is in the 
fulfilment of the following requirements through the above described 
stages and all variations/modifications thereof: 

10 1) Sample preparation 

(definition of the Afi ch for a given material (sample)) 
2) Experimental arrangement 

(definition of the sample geometry within the given 
electrical measurement circuit) 
15 3) Method of measurement 

(measurement of the electrical response) 
4) Method of analysis 

ch 

(calculation of the electrical material parameters Lfi , 
e 9 N and a (here also all possible combinations of these) 
20 by using various analytic expressions for parts or the 

whole of the measured Z (w) function or by using the 
numerical solutions of the governing equations with the 
electrical material parameters as free fitting parameters. 

25 In what follows the individual stages 1) to 4) will be described in 
detail using figures 1 to 5. The new method for electrical 
characterization of semiconducting/insulating materials will be 
demonstrated on a couple of examples shown in figures 6 to 8. 

30 Finally Table 2 illustrates the quantitative results for the values 

rh 

of A/i , c 9 N and a on a number of monocrystalline silicon samples 
of both n and p-type with various degrees of doping, obtained by the 
use of the present invention and compared with the values for these 
parameters obtained by other methods on the same samples and from 
35 Physics tables (4) and the literature. 

stage i - Sample preparation 

(definition of Lfi ch = M° h B - M Ch M+s of the material) 
In order to achieve the first requirement (establishment of the 
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conditionUp^Uo for the material) reference is made to Fig. 1 for 
a sample sandwiched between two metal electrodes in a preferred 
geometric configuration depicted in Fig. 2a (a simple capacitor 
configuration). The problem in this case is quasi one-dimensional 
5 (area A p of the sample is constant and its linear dimension Jh is 
large compared with the length of the sample L). 

Referring to Fig. 1 the sample of the finite length L extends from 
x=0 to x»L, x=0 and x=L planes being the true surfaces of the 
10 sample. The metal electrodes are deposited Ax --> 0 (or applied Ax > 
0) on both these surfaces (capacitive coupling of the sample to the 
external field) and when an external electrical voltage is applied 
between these two electrodes an external electrical field is created 
within the volume of the sample. 

15 

With no external electrical voltage applied the average energy of 
the mobile electrical charges well within the bulk of the sample is 

eh 

characterized by the thermodynamical quantity fig - chemical 
potential of the bulk. The average energy of the electrical charges 
20 at the surface of the sample (and/or of sample surface-metal 
electrode interface) is characterized by the sample surface-metal 
electrode interface chemical potential M^^s* 

Considering the n-type material first there will be an outflow 
25 (inflow) of the mobile electrical charges from (into) the sample to 
(from) sample surface-metal electrode Interface region depending on 
whether the chemical potential /ig Ch of the bulk of the sample lies 
above (below) the value of the chemical potential of the sample 
surface-metal electrode interface region M 0 * 1 ^* 

30 

This flow of the mobile electrical charges from the sample interior 
(sample exterior) creates an internal electrical field in the 
regions near the surfaces of the sample that eventually stops any 
further flow. In this way a dynamical thermodynamical equilibrium is 
35 established that is characterized by the existence of the 

depletion regions (outflow of the mobile electrical charges from the 
interior of the sample into the sample surface-metal electrode 
interface region or 

accumulation regions (inflow of the mobile electrical charges from 
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the sample surface-metal electrode interface region into the 
interior of the sample)* 

The requirement of the Stage 1 of the present invention (|Am°^I> 0) 
5 is simply the requirement for establishing (securing the existence 
of) the depletion regions in the material by securing that for the 
n-type material fi c ^ D > p c ^ u > c (the same arguments apply for the 

rh rh 

p-type material but here the condition is n M+s > 0 

10 If the condition p ch B > l» ch || fS (n-type material) or M Ch M+s > /* ch B 

(p-type material) are not satisfied then an appropriate treatment of 

the sample surface might be necessary. This applies in cases where 

M^N+j is determined mainly by the surface states of the material. 

Alternatively an adequate metal electrode - with the right value of 

15 the work function - has to be chosen. This applies in cases where 
ch 

N+S is determined mainly by the work function of the metal 
electrode. In general a combination of both processes might be 
required. 



20 Both of these preparation processes or steps move the value of 
ch 

fi M+s up or down on the energy scale (see Fig. 1) so that the 
condition Ap 0 * 1 > 0 can be eventually obtained for n-type and p-type 
materials (A/i ch < 0). 

25 Although a specific description of the Stage 1 has just been given, 

it is clear to anyone experienced in the field that other ways of 

ch 

securing the condition A/r > 0 might be possible. 



The Stage 1 of the present invention therefore does not relate to 

30 any particular process or sequence of processes how to achieve the 
ch 

condition Am 0 but it is simply embodied in the condition itself 

ch 

- establishment of the A/r 0 in the material, (for n-type and 
p-type material, respectively.) 



35 Stqqe 2 -Experimental arrangement 

(definition of the sample geometry, of the sample 
environment and of the electrical measurement circuit) 



Experimental arrangement that constitutes the Stage 2 of the present 
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Invention is defined in three steps. The first is sample - metal 
electrodes configuration, the second is definition of the sample 
geometry in a given sample - metal electrode configuration and the 
third step defines the sample environment together with electrical 
5 connections to the external electrical measurement circuit. These 
steps will now be described in more detail. 

Sample - metal electrodes configuration 

defines the way the external electrical field is coupled to the 
10 piece of the material under investigation. This can be done in a 
variety of ways in principle but in the preferred embodiment 
(preferred experimental arrangement) a most simple capacitive 
coupling is preferred where the material and the electrodes form a 
simple capacitor. 

15 

In this configuration, known also as "sandwich configuration", the 
electrical field profile within the sample is simple, making the 
definition of the sample geometry in the electrical measurement 
circuit easy and precise. 

20 

Two basic preferred configurations are considered. They are both 
shown in Figures 2a and 2b. 

In the "research" configuration (Fig. 2a) the appropriate metal 
25 electrodes are deposited directly onto both surfaces of the material 
(Ax --> 0 in Fig. 1} minimizing the "edge" effects of this simple 
capacitor (no corrections for the "edge" effects in the calculation 
of the sample geometry is needed). 

30 in the "on-line quality control" configuration (Fig. 2b) the major 
requirement is that the electrodes can be applied at a particular 
measurement point Pj (Xj,yj) of the semiconductor wafer surface and 
after the measurement of the electrical response is completed at 
this point the electrodes can be re-applied at a new point ?2^ x 2' y 2^ 

35 with a new measurement of the electrical response to follow. The 
process is repeated for all the points P^(x^,y^) of the wafer 
surface. The effective electrode area (including the "edge" effects) 
is required to be the same through-out the whole of the mapping 
measurement process (measurement of the electrical response across 



WO 94/14078 



PCT/DK93/00410 



14 

the whole area of the wafer). 

Another possible sample - electrode configuration 1s shown in Fig. 
2c- This type of configuration - an "open capacitor" configuration 
5 (also known as "planar" configuration) might be relevant for the 
measurements on thin films and/or semiconducting/insulating bulk 
materials (here also wafers) where a simple capacitor configuration 
is not possible. The preferred shape of the electrodes in this case 
will be rectangular with well defined gap between them. 

10 

In all the configurations described above the geometrical shape of 
the electrodes 1s not Important but must be known together with the 
corresponding external electrical field profile. 

15 Sample geometry 

In a given sample - electrode configuration, either a preferred 
embodiment - Fig. 2a and 2b or any modification thereof such as 
shown for example in Fig. 2c, the next step in Stage 2 of the 

20 present Invention consists of the determination of the geometrical 
dimensions of the system sample + metal electrodes that are relevant 
for the definition of the given configuration as a electrical 
element within the electrical measurement circuit. The quantities in 
question are the areas and the lengths (thicknesses) of the 

25 individual elements comprising the sample - electrode configuration* 

Referring to Fig. 2d (the "research" configuration), the elements of 
the configuration are the top metal electrode (of area A f and length 
tj), the piece of the material under investigation (of area A r and 

30 length L) and the bottom metal electrode (of area A f and length t 2 ). 
If the electrical resistance of the metal electrodes is low (low 
compared to the resistance of the piece of the material under 
investigation), the effect of the metal electrode elements on the 
electrical response can be (safely) neglected and the geometrical 

35 dimensions of the sample - electrode configuration is that of the 
sample Itself (area A f and the length L). 

In the "on-line quality control" configuration shown in Fig. 2e the 
same arguments apply and the relevant elements (with their 
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geometrical dimensions) are the contacting conductive medium (of 
area Aj and the length A x; in the preferred embodiment of the 
present invention these quantities are the same for both electrical 
contacts), the piece of the material under investigation (of the 
5 effective area A ? (inclusive the "edge" effects) and the length L) 
and the top and bottom metal electrodes (both of the area A3 and the 
length t). 

All these quantities have to be measured with a sufficient precision 
10 prior to calculation of the four electrical material parameters 
Aji ch , €, N and n since the relative errors in their measurement 
determines the overall accuracy with which the numerical values of 
the above mentioned four electrical material parameters can be 
obtained using the present invention. 

15 

Sample environment and electrical contacts 

Two basic sample environments are envisaged in the preferred 
embodiment of the present invention. 

20 

In the "research" configuration the sample (F1g. 2d) is placed in 
the cry os tat, so that the dependence of the four electrical material 
parameters A/i ch , c, N and m can be studied as function of 
temperature, pressure, d.c. electrical bias, optical (particle) 
25 excitation of, the electrical charges within the material under 
investigation etc. 

In the "on-line quality control" configuration, the sample (a 
semiconductor wafer - see Fig. 2e) 1s placed preferably In a 

30 automatic wafer exchange mechanism that introduces one wafer at a 
time into a dark (optically) measurement compartment, where x,y 
mapping of the four electrical material parameters takes place at 
room temperature, the whole assembly making the large throughput 
(number of wafers processed in a unit of time) possible. Other 

35 sample environments are evidently also possible and this step of the 
Stage 2 of the present invention does not therefore relate to any 
particular sample environment but only demands that the sample 
environment be well known (temperature of the sample, definition of 
the "optically" dark measurement compartment etc.). 
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Once placed in its appropriate environment the sample has to be 
connected via electrical leads to the external measurement circuit. 
In the preferred embodiment of the present invention this is 
achieved best by connecting four screened coaxial lines (two and two 
5 connected together at the sample - see Figs. 2a-2c) to each of the 
sample surfaces (each of the sample electrodes), minimizing in this 
way the unwanted effects of the external leads (the self-induction 
and the stray capacitances of the electrical leads) on the measured 
electrical response of the sample* Two of the leads (one from each 
10 side of the sample) are used to measure the electrical voltage 
across the sample while the other two leads are used to measure the 
electrical current through the sample. The other ends of the four 
coaxial lines are then connected to the external electrical 
measurement circuit. 

15 

Other possible ways of connecting the sample to the external 
electrical measurement circuit exist evidently also. 

This step of the Stage 2 of the present invention (electrical 
20 contacts to the sample) does not relate to any particular way of 
connecting the sample to the external electrical measurement circuit 
but resides (is embodied) in the requirement that the effect of the 
external electrical leads to the sample in the alternatives 
(modifications) to the preferred embodiment (Fig. 2a and Fig. 2b) 
25 must be known and properly accounted for in the final analysis of 
the electrical response of the whole system, the system being 
defined as consisting of the sample, the electrodes and the 
electrical leads. 

30 Stage 3 - Method of measurement 

(definition of the measurement method) 

There exists a number of possible ways of measuring the electrical 
response of a system (sample, electrical metal contacts and external 
35 electrical leads) involving either real time and/or frequency 
domain. All these experimental methods involve a measurement of the 
two relevant quantities - incut signal (be that electrical charge, 
electrical voltage or electrical current) and the 
- output signal (complementary quantity to the input), the ratio 
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between the output and the Input defining the appropriate electrical 
response function. Two forms of the input signal are used most 
commonly. These are the "step" or "sinusoidal" input. The respective 
electrical response functions of these two input forms are related 
5 to each other through Laplace transform (linear response regime). 

A typical electrical measurement arrangement is shown schematically 
in Fig. 3. In the preferred embodiment of the present invention a 
sinusoidal voltage from the voltage source is applied across the 

10 sample (voltage input) and its amplitude is measured by a "vector" 
voltmeter. The amplitude and the phase shift (relative to the 
voltage input) of the output current (electrical current measured in 
the external electrical circuit) 1s measured by a "vector" ammeter 
and the ratio of the in-phase and out-of -phase component of the 

15 current output to the voltage input then defines the complex 
electrical admittance Y (w) of the system under investigation 
(system being defined as sample + electrodes + external electrical 
leads). 

20 Another quantity that characterizes the electrical response of the 
system under investigation is the complex electrical impedance of 
the system defined as Z («) = 1/Y (w). 

The direct experimental determination of the Z (w) is accomplished 
25 in practice most conveniently by various types of low frequency 

(w < 1.10 10 Hz) impedance a.c. bridges consisting of a voltage 
source with a variable frequency f (a = 2Hf), preferably with a 
range of 1 mHz through 10 GHz, "vector" voltmeter, "vector" ammeter 
and some additional circuitry to yield the function Z (t>) at various 
30 frequencies of the input sinusoidal voltage. All these components 
form elements of a single measurement unit shown schematically in 
Fig. 2a and Fig. 2b. 

The Stage 3 of the present invention does not depend on any 
35 particular way of the measurement of the electrical response of the 
system under investigation but relates to a process consisting of 
the measurement of the electrical response of the system and its 
subsequent representation in the frequency domain by the complex 
electrical impedance (admittance) Z (w) (Y (w)). This representation 
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is not (strictly) necessary and other representations of the 
electrical response of the system under investigation, such as for 
example the complex dielectric constant T (0), the complex 
electrical conductivity a (u), the complex electrical modulus M (w) 
5 and /or various variations of the electrical response functions in 
the real time domain, are all of course possible. The Z (w) (Y (0)) 
representation of the electrical response of the system is preferred 
for the reasons of clarity of the physical arguments behind the 
analysis of the electrical response data and because the numerical 
10 values of these functions at various frequencies have a very 
definite physical meaning in relation to the four electrical 
material parametei 
such an analysis. 



25 



rh 

material parameters A/i , £, N and 0 that are to be extracted from 



15 Within the necessary experimental range of frequencies f (1 mHz - 10 
GHz) will these numerical values fall within the following typical 
ranges for most of studied semiconductors/insulators: 

14 

the real part of impedance in the range from 0.001 ohm to 10 ohm. 

-14 

the imaginary part of admittance divided by u from 10 Coulomb to 
20 10" Z Coulomb. 

Stage 4 - Method of analysis 

(calculation of the electrical material parameters 
Ln , c, N and m of a material under investigation). 



Introduction 



The stage 4 of the present invention is the crucial step. It is at 
this stage that the known physical theories and their mathematical 

30 formulation are applied to a specific technical problem, namely the 
problem of the electrical characterization of semiconducting and/or 
insulating materials in general ("research" configuration of the 
present invention) and in a more specific application, the problem 
or the electrical characterization of semiconductor wafers ("on-line 

35 quality control" configuration of the present invention). 

The discovery aspect of the present invention can be summarized in 
two major points: 
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Firstly the discovery of a direct link (that can be formulated 
simply and clearly through mathematical equations and/or 
expressions) between the full theoretical analysis of the electrical 
response of a material, including the depletion regions and material 
5 surface-metal electrode interface regions and the experimentally 
observable electrical response function Z (o>) of the material (where 
the function Z (w) describes the dependence of the complex 
electrical impedance of the material on the frequency of the applied 
external voltage input). 

10 

Secondly the discovery of a possibility of extracting the four 

ch 

electrical material parameters A/t , c, N and fi of the material 
(and/or their combinations) when the depletion and material 
surface-metal electrode interface regions are included in such a 
15 theoretical analysis. 

The invention itself is then embodied in the application of this 

discovery aspect to the technical problem of characterizing the 

electrical properties of semiconducting/insulating materials and 

20 semiconductor wafers through the determination of the four 

ch 

electrical material parameters A/T \ c 9 N and fi using the above 
mentioned method of analysis of the experimentally observed function 
1 (ft)). One of the claims of the present invention will be therefore 
the use of the present invention to determine one, two, three, four 
25 and/or any combination of the four, electrical material parameters to 
characterize electrically a piece of a material. 

Because of the importance of the Stage 4 of the present invention a 
more detailed description of this stage will now be given. 

30 

Firstly, the definition of the problem will be presented. 

Secondly, the description of the essential parts of the theory 
behind the full theoretical analysis of the electrical response of a 
35 real material (material with depletion regions, surface and material 
surface-metal electrode interface regions) will follow. 

Thirdly, the various possible approximations and the corresponding 
solutions will be described and the link between the governing 
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equations (and the corresponding solutions thereof) and the 
experimentally observable response function 1 (w) will be pointed 
out. 

5 Fourthly the various methods of calculating the numerical values of 
the four electrical material parameters A/i c *\ £, N and /x will be 
described and the validity of the method of analysis of the present 
invention will be demonstrated on a set of experimental results 
obtained by the use of the present invention on a number of 
10 monocrystalline and polycrystalline silicon samples of both types 
(n-type and p-type) and with varying degrees of doping. 

These results (the obtained numerical values of the four electrical 

ch 

material parameters A/r , £, N and p) will be then contrasted in the 
15 final Table 2 with the results for the same four electrical material 
parameters obtained either on the same samples by the other 
experimental methods or by comparing the obtained values with those 
in the literature and/or in Physical Tables of material constants 
(4). 

20 

Definition of the problem. 

Referring to Figs. 2a, 2b and 3 the measured electrical response 

e*+ 

1 (a>) is due to the electrical response of the sample (bulk region, 
25 depletion regions and the surface-metal electrode Interface regi- 
ons), the electrical response of the contacting, conductive medium 
(Fig. 2b), the electrical response of the metal electrodes (Figs. 2a 
and 2b), the electrical response of the electrical contacts (Figs. 
2a and 2b) and the electrical response of the external electrical 
30 leads (Figs. 2a and 2b) connecting the sample configuration within a 
given sample environment to the external measurement circuit. 

In order to obtain the electrical response function (w) of the 
sample alone, it is therefore necessary to subtract from the total 
35 electrical response function lj Q j (<u) of the total system the 
electrical response from all the unwanted elements of the system 
(contacting, conductive medium (if present), metal electrodes, 
electrical contacts and the external electrical leads. 
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This subtraction is in principle easy since the elements of the 

system are connected in series so that the total electrical response 

function Zj oT (a>) of the system is the sum of the complex electrical 

impedances of the individual elements of the system. Provided that 
_ #v 

5 the electrical impedances Z i (w) of the individual elements of the 
system are known (this can be rather difficult to achieve in 
practice), this subtraction can be performed. 

In the preferred embodiment of the present invention, the effect of 
10 i (o>) of the elements of the system on the resulting total 
electrical response function Z ToT («) 1s minimized, so that to a 
very good degree of precision (typically better than 1% accuracy) 
these effects can be neglected and the measured total electrical 
response function Zj Q y (w) of the system represents the true 
15 electrical response of the sample Z s (w). 

Having extracted the electrical response function zl (o>) of the 
sample from the experimentally observed lj Q j (w) (if modifications 
to the preferred embodiment of the present invention are used), the 
20 next step is to subject this function to the theoretical analysis 
embodied in the present invention. 

The application of this analysis then results in the possibility of 
the determination of the numerical values for the four electrical 
25 material parameters A/i ch , £, N and fi (and/or their combinations such 
as for example the specific resistivity of the bulk region of the 
material) that characterize the electrical properties of the 
material (semiconductor, insulator, semoconductor wafer etc*) 

30 The results of the present analysis can be summarized very 
schematically as follows: 

By including depletion regions and the surface (and/or material 
surface-metal electrode interface) regions of the material under 
35 investigation in the full analysis of the electrical response of the 
sample (rather than to try to avoid it at all costs, as has been the 
case up till now), the analysis method of the present invention 
shows that the form or the values of Z $ (g>) in various frequency 
regions and/or at various frequencies are related directly to the 
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ch 

four electrical material parameters Ap , e, N and of the 
material . 

This means for example that under appropriate experimental 
5 conditions the value of the real part of the Z $ (w) at high 
frequencies is related directly to the bulk resistivity p of the 
material (p = 1 / N . e . m)> while the imaginary part of the Y $ (w) 
divided by 0 is related directly to the bulk dielectric constant e. 

10 Similarly, at low frequencies both the real part of Z c (*>) and the 

imaginary part of Y r (a>) divided by cj are directly related to the 
ch 

parameters A/i and N. 

Furthermore, by looking at the form of Z s (w) (and/or Y $ (w)) in the 
15 appropriate frequency regions it is possible to establish the type 
of the majority charge carriers and whether the electrical transport 
within the material is due to only one kind of mobile electrical 
charges (electrons or holes) or whether both mobile electrical 
charges (electrons and holes) contribute significantly. 

20 

Finally, it is possible for the first time to study quantitatively 
the bulk region, the depletion regions and the surfaces (and/or 
material surface-metal electrode interface regions) of the material 
under Investigation Independently and separately by measuring 
25 (studying) the electrical response function Z (id) in various 
frequency regions. 

Description of the theory behind the analysis method 
(the elements off) 

30 

a) Electrical response in solids at classical frequencies 

A new analysis of the electrical response of a macroscopic solid 
state system to an arbitrary electrical voltage input indicates that 
35 with a small sinusoidal voltage input within the classical frequency 
range (w < 10 10 Hz), the electrical response In a solid state system 
mirrors first of all the static, spatial distribution of the mobile 
charge carrier density through the local electrical conductivity 
relaxation time T ff (x) = c(x)/a D c (x) (c (x) is the local dielectric 
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constant, a D c (x) is the local d.c. electrical conductivity) (see 
Figs. 4b and 4c). 

This enables any system, inclusive the depletion regions and the 
5 electrical contact regions of a material to be modelled by a simple, 
passive R,C electrical network, where each of the network elements 
has a very direct physical meaning. 



Furthermore, since the electrical conductivity relaxation time is 
10 more or less the only relevant characteristic time of the problem, 
the regions of the sample with different r Q will respond at 
different characteristic frequencies that are directly related to 
r a . In this way it is therefore possible to study for example the 
bulk, the depletion and the electrical contact regions separately by 
15 studying the electrical response of the sample (system) in different 
frequency regions. 

b) Mathematical formulation 



The macroscopic electrical response of a medium (characterized by a 
dielectric constant t and the electrical mobility p) at classical 
frequencies is described fully by classical electrodynamics (two of 
Maxwell equations), the constitutive equation defining the total 
local electrical current, the initial condition for the total, 
electrical charge density p(x,t) spatial distribution and by the 
boundary condition defining the charge transport across the 
boundaries. 

Within the classical range of frequencies (a> < 10 10 Hz) the 
30 dielectric constant c of the material and the electrical mobility m 
of the mobile electrical charges within the (bulk of) the material 
can be both considered as time independent (possible slow polarisa- 
tions neglected in the first approximation). 

35 In one dimension (the appropriate experimental arrangement) and when 
both types of the mobile charge carriers (electrons at E c and holes 
at E y - see Fig. 1) contribute to the transport, the defining 
equations can be re-cast into a set of coupled, non-1 1near parabolic 
equations for each type of the mobile charge carrier particle 



20 



25 
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density: 



5 =. fc .m*±iL^ (Ec , x , a . KA x. a .*.t± ;; *.a 

M ax 2 • * 



10 



15 



20 



~T3T* * - R h 



d KuxtEj . x.C) 

Tt = -| (E.R.), - (C.J?.), 



= -[(£.«.), - <C.tf.),i £/ 



(1) 



where 

25 

div t) = Q (x * c) 
c 



Elx.O = E rwr (x.C) ♦ ^(t) < 2 » 

35 
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5 



The constitutive equation defining the total local elec 
trical current has the form 



jix, t) = |e| • [u e - n 0 (E c ,x,t) + \i h - n h (E v ,x,t)] • E(x, t) - 



jc-T f dn e (E c ,x,t) BnAE v ,x, t) 1 



10 



and the boundary condition for e(x,t) is defined through the 
relation 



Here p c and m are the electrical mobilities of the elec- 
trons at E c and of the holes, at respectively, c is the 
20 dielectric constant of the material, k is Boltzmann constant, 
T is the absolute temperature and |e| is the elemental elec- 
trical charge* 

The quantities n e (E c ,x,t), nJE^x^), E (x,t) and Q(x,t) 

are the particle number densities of the respective mobile 
charge carriers, the local electrical field and the local, 
25 total charge density respectively- N t mUK S(E t ,x, t ) and 

Rj ,u>c"(E J ,x, t) are the local particle number densities of char- 
ged, localised energy levels E t and within the energy band 
gap E g - 

(E,R. )i , (C,R, ) 4 , (E.R. )j and (C.R. )j are their respec- 
tive emission and capture rates and finally G e , R e and G h , R h 
are the generation and recombination rates for mobile 
electrons at E c and holes at E^ respectively. 



15 




(4) 
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c) Results - Dynamical solution 

A set of the governing equations (1) and (2) together with the 
appropriate initial and boundary conditions for n e (E c ,x,t), 
5 n.(E v ,x,t) (mobile electron and hole particle densities), 

N 1 L0C (E r x,t) and N "j,L0C tE j ,x,t) ( localized Positive and negative 
charge particle number densities) then determine the space-time 
evolution of the total charge density p(x,t) and therefore also the 
electrical current flowing 1n the external circuit in a response to 
10 a given applied voltage input, thereby defining the electrical 
complex impedance of the sample Z $ (w). 

d) Results - Long time, static limit 

15 For times t » 0 (t--> *) and with no external applied field, the 
set of the governing equations leads to the formation of the 
(majority) mobile charge carrier depleted regions near the surfaces 
if the boundary condition reflects the finite difference in the 
chemical potential across the boundary at time t « 0. In this case 

20 the space-time evolution of the total charge density />(x,t) will 
approach the equilibrium, time independent distribution as the time 
goes to "infinity" t Q — > *. This distribution 1s also a solution of 
the static Poisson equation to which the problem reduces in this 
long time, static limit. 

25 

e) Results - Small signal approximation and the static R.C network 
solution to the problem of the elec trical response 

One of the major results of the present dynamical analysis of the 
30 electrical response in sol Ids at classical frequencies is concerned 
with a small input signal approximation. 

When the external applied voltage is sufficiently small, the 
equilibrium, spatially non-homogeneous distribution of the mobile 
35 charge carrier particle densities will not be disturbed by the 
applied external electrical field and it can be shown that the 
electrical response of the entire system (sample) under these 
conditions is identical to a response from simple, parallel R,C 
electrical elements, connected in series. 
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Contrary to the usual passive R,C network models of various 
junctions, interfaces and semiconductor - insulator - metal 
structures, the electrical elements in the present static R,C 
network have a very direct physical meaning and are all 
5 interrelated* 

According to the result of the present analysis, the sample is 
simply divided into a number of volume elements (see Figs. 4a, 4b and 
4c), the actual number depending on the required precision with 

10 which the electrical response Is required. Each volume element AV,- 
is characterized by Its electrical resistance Rj (in-phase component 
of the response; usually dissipation of the external field energy 
through finite electrical mobility /i of the material) and by its 
geometrical capacitance Cj (out-of phase component of the response; 

15 fast polarisation through finite dielectric constant c of the 
material). 

The solution of the Poisson equation in the long time, static limit 
determines the spatial distribution of the electrical mobile charge 
20 density within the sample and therefore also the electrical 
resistance of each of the volume elements AV (Xj). The geometrical 
dimensions of each volume element then determine its characteristic 
geometrical capacitance. 



25 By defining R's and C's for all the volume elements comprising the 
sample (bulk region), the depletion regions and surface-metal 
electrode interface regions) one has determined also the total 
electrical complex impedance of the sample Z $ (w) which is simply 
the sum of the impedance of the individual volume elements. 

30 

As can be seen from Fig. 4c each of the volume elements has got its 
own characteristic electrical conductivity relaxation time Tj and 
the "visibility" of the depletion regions, bulk region and the 
contact regions of the sample in the experimentally observed Z $ (w) 
35 becomes now obvious. Depending on the frequency of the applied 
external electrical field (applied external voltage input), the 
different parts of the sample will be (or will not be) able to 
respond. 
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In the long time, static approximation therefore, the measurement of 
the electrical response (Z $ (0)) at different frequencies is 
directly related to the equilibrium spatial distribution of the 
mobile electrical charge particle density in the different parts of 
5 the sample. 

Identification of the link (correlation) between the theoretical 
analysis of the ele ctrical response and the experientallv observed 
electrical response function Z $ (w). 

10 

It is the discovery of a direct correlation between the space-time 
evolution of the total charge density p(x,t) within the sample and 
the experimentally observable electrical impedance of the sample 
Z s {*>) together with the solutions for p(x,t) (either a dynamical 

15 solution obtained by solving the whole set of governing equations 
numerically or a static solution obtained by solving the static, 
time independent Poisson equation with corresponding reticulation of 
the sample into the sum of the volume elements, each characterized 
by its own impedance), that constitutes the discovery aspect in the 

20 theoretical analysis embodied in the present invention. 

The use of this discovery aspect to characterize the electrical 
properties of semiconducting/insulating materials and semiconducting 
wafers through the extraction of the four electrical material 
25 parameters Ap ch , c, N and n (either some of them, all of them or 
some combination of them) using this method of analysis (together 
with the other stages described) then constitutes the present 
invention. 

ch 

30 Calculation of the four electrical material parameters AiT 1 . £. N 
and u of the material under investigation. 

Depending on the particular material under investigation, the 
calculation of the electrical material parameters from the 
35 experimentally observed electrical response function (complex 
electrical Impedance) of the sample Z $ (<■>) can proceed (can be 
performed) at three possible levels described in the foregoing 
paragraphs. These levels will now be described in more detail in 
order of increasing complexity. 
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Level 1 - "Hat" model approximation. 

It turns out that within the long time limit static approximation a 
further simplification can be employed in the static R,C network 
5 analysis of the electrical response of the sample that makes the 
determination of the four electrical material parameters Aji ch , c, N 
and n even more direct. If the (characteristic) frequency regions 
corresponding to the response from the material surface-metal 
electrode interface regions, from the depletion regions and from the 

10 bulk region of the sample are all well separated, it is possible to 
lump all the volume elements of the bulk region together and do the 
same for the depletion regions and the sample surface-metal 
electrode interfaces. The electrical response of the sample is then 
represented by a very simple static R,C network consisting now of 

15 only three R,C elements representing in turn the electrical response 
from the sample surface-metal electrode interfaces (time constant 
r H+s ), the electrical response from the depletion regions (time 
constant 7 Q ) and the electrical response from the bulk of the sample 
(time constant r B ). This representation (named here as "Hat" model 

20 approximation) of the electrical response of the sample shown in 
Fig. 4d is sufficiently precise only in case when the characteristic 
time constants of the respective regions of the sample are well 
separated (^5 > > 7g). 

25 The "Hat" model reticulation of the sample shown in Fig. 4d enables 
one to write down analytic expressions for the four quantities Rp, 
C D , R B and C g that can be identified in the "raw" experimental data 
of the measured Z s (u) of the sample. 

30 The expected frequency dependence of Z $ (a>) for a sample represented 
by the equivalent electrical circuit of Fig. 4d is shown in Fig. 5. 
For clarity purpose the real part of (w) (dimension Ohm) and 
imaginary part of y! (g>) / 0 (dimension Farad) are displayed. The 
flat parts of the curves denoted by symbols R^y R^ c (g) and c (5) 

35 are directly related to R Q , Rg, C D and C g through algebraic rela- 
tions 
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R m = 2 -«„(D 



10 



R W ~ R i 



(1) 



(T) 



'(2) 



1 - c D (D 



(5) 



(3) 



= C fl (T) 



15 



which in turn are defined in terras of the four electrical material 
ch 

parameters Ap , c, N and ft of the material under investigation. 

20 For the case of a typical semi conductor/ insulator with a band gap E g 
and with one-particle effective density of states N c (conduction 
band) and N y (valence band) the defining expressions are 



25 



30 



35 



R 1 • L " 2 ' Wd{T> 

B ~ |e|-n(D-w A r 
C B = c 0 -c c -- *'„ IM (6) 



L - 2 W_(D 



e »** . 2-A t i rfl (D 



arc tan 
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where 



u = e 



10 



vr n = 



c = 



c = 



= 1 2 *si 



ch 



and a = Etf+An** 



and a = JSJ s p -A|i cb (p- type) 



15 



P 2 = i^ 



The quantities appearing in the analytical expressions 
2q have the following meaning : 

P (T) - electrical mobility of the mobile, majority $1. 

charge carriers 
P C (T) - electrical mobility of the electrons at E c 

P h (T) - electrical mobility of the holes at E v 

„ N - mobile, majority el. charge carriers particle 

" number 

density 

L - total length of the sample 

W 0 (T) - depletion width of the sample 

K - total, active area of the sample 

£ 0 .e r =:£ - dielectric constant of the sample (silicon) 

E fl - energy gap of the material 

30 Ap^T) - difference in the chemical potential 

N c ,Nv - one-particle effective densities of states in 

the conduction (Sc) and valence (Ey) bands 
respectively 
\e\ - elemental electrical charge 

k - Boltzraann constant 

T - absolute temperature 

35 
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Considering now the four electrical material parameters A/* c \ c, N 
and ix as four unknown constants defined by four independent 
equations, it is clear that a unique solution exists. It can be 
found iteratively by first assuming some definite values (starting 
5 values) for A/i ch , c, N and then calculating the quantities 

R (D* C^and comparing the result with the observed values of 

R R C 0 and C ® ™ e next step in tlle iterat * on procedure 
consists in changing the starting values for A/x ch , c, N, /t and 

repeating the whole process until the required precision (agreement 

10 between the calculated and measured values of C^and 

is obtained. The fitting procedure is so simple that it can be done 

by hand and does not require computer. In these cases where the 

ratio of the electrical mobilities P ? - Mj/Mg 1s also considered as 

an unknown parameter, a combination of the measurements on both 

15 n-type and p-type samples 1s needed in order to obtain a unique 

sol ution. 



Level 2 - Static approximation using the solution of the 
20 Poisson equation 



When a fit to the experimentally observed electrical response 
function Z $ («) of the sample over the entire frequency domain 

25 measured in the experiment is required for higher precision or when 
the individual levels R^ R q C ^ and C q can not be easily 
identified in the raw experimental data Z (w), the static R,C 
network analysis of the electrical response of the sample involves 
the reticulation of the sample as shown in Fig. 4e. Here the 

30 depletion region volume elements have not been lumped together but 
rather the resistances of the individual volume elements within the 
depletion region have been determined through the solution of the 
Poisson equation for the total local charge density p (x,t > ») 
assuming simple parabolic bands (parabolic band approximation - see 

35 Fig. 1). The electrical impedance of the depletion region can then 
be calculated analytically leading to the expression (both electrons 
and holes contributing to the electrical transport) for Z Q (w) of 
the form 
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10 



Eg 



arc tan 
k-T 



15 t x = t 0 • e 



2-k-T 



20 



P' = P 0 •( e * r - l ) 
t = P 0 •( + i ) 



2-A|i 



ch 



( 1 + i • ot) / 



*0 = 



2-|e| • JN C - N v - p 0 > H (7) 



Po - 



a 



-.5. 



where 



25 



a = 



30 
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In order to extract the four electrical material parameters Aji Ch , c, 
N and fi of the material under Investigation, the same iterative 
procedure as before is used only now the fit can be best performed 
by a computer using the four electrical material parameters as free 
5 fitting parameters. The computer is used mainly for the purpose of 
time economy, since in principle the calculation can be done by hand 
due to the fact that the solution for the Zg (a) is analytical also 
1n this case. If the analytical expressions for the interface 
element (material surface-metal electrode interface region) are not 

10 known, the fit to the experimentally observed Z Q (0) 1s performed 
using and C M+S as two extra phenomenological free fitting 

parameters. In case of silicon though the analytical expressions for 
the interface elements R^, C M+S do exist as will be demonstrated 
in the experimental verification of the present invention. The 

15 effect of and beyond the frequency 1 / r M+$ (Fig. 5) Is 
small, the form of the response function Z $ (w) in this frequency 
region being determined mainly by Z* D (w) and the known bulk element 
Z B (u) (elements R g and C g ). 

20 Level 3 - Dynamical solution 

When a complete analysis of the electrical response of the material 
1s required either for reasons of highest precision or when the form 
of the Z $ (to) is more complex and the detailed analysis of the 

25 boundaries of the sample (interface element R^^ and C M+$ ) is 
needed, the full dynamical solution of the problem has to be 
undertaken numerically. This situation corresponds to Fig. 4a and no 
reticulation of the sample into individual volume elements is 
performed (reticulation is not relevant in this case since the 

30 reticulation of the sample and the static R,C network analysis of 
the electrical response of the sample is the approximative solution 
representation of the full dynamical solution of the problem). 

Besides the governing (transport) equations for the local charge 
35 density p(x,t) space-time evolution describing the electrical 
response of the sample to the external voltage input, 1t is 
necessary to define (to model) the initial condition for the p (x,t) 
and the appropriate boundary conditions (transport of the electrical 
charge across the boundaries - surface of the sample). 
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Once the Initial and the boundary conditions are defined the set of 

the governing equations for the space-time evolution of the total, 

local charge density p (x,t) is solved numerically with some 

ch 

starting values of the four electrical material parameters AfT , c, 
5 N and p. The resulting calculated electrical response function 

Z CAL of the sam P 1e is then compared with the experimentally 
observed complex electrical Impedance Z^ xp («) and the whole process 
undergoes iterations until the required agreement between Zq^ (g>) 
and Z EXp [w) is obtained, leading to the final values of the four 
10 electrical material parameters Am C **» c» N and n characterizing the 
electrical properties of the material under Investigation. 

At this level of solution the computer 1s essential since the 

+++ 

problem of calculating (a) can be done only numerically. 



15 



Verification of the validity of the present invention 
(Experimental examples) 



The first set of experimental tests has been performed with pure, 
20 monocrystalline silicon samples of both types (n-type and p-type) 
with varying degree of doping using either evaporated thick gold 
film or a conductive silver paint as metal electrode material. 



Typical raw experimental results of the electrical response in these 
25 samples are shown in Figs. 6a and 6b where log (Re (Z $ (w))} and 

log (Im (Y s (w) / w)} are plotted versus log of the measurement 
frequency for a number of different temperatures. 

The verification of the validity of the present invention has been 
30 proved In two different ways. The first 1s demonstrated in the Table 
2 where the numerical values of the four electrical material para- 

ch 

meters A/i , t f N and p are compared either with the literature 
values (A/x ch > c and n) or with the values obtained by using 4PDCR 
and FTIR/FTIPL experimental methods on the same samples. 

35 

/A high degree of agreement has been obtained (to within 5%) limited 
only by the precision of the determination of the geometrical 
dimensions of the sample. 
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As is apparent from the Table 2 the values of the four electrical 
material parameters determined by the "New" EIS method (present 
invention) are more precise than the values for these parameters 
obtained by the other methods or obtained from the literature. 

5 

The second way of verifying (testing) the validity of the present 
invention was the use of the temperature variations of the 
experimentally observed Z s (*>) of the samples (temperature being 
used as a known parameter). 

10 

The four electrical material parameters have been extracted from the 
measured Z $ (w) at a given temperature by fitting these four para- 
meters at this single temperature. The form and the numerical values 
of the Z Q (u, T) were then calculated at other chosen temperatures 

15 (the temperature dependence of the Z<* (u, T) is known - it is also 
one of the results of the theoretical analysis of the electrical 
response embodied in the present invention). The calculated values 
for Z s (w, T) were then compared with the observed Z $ (w, T) at 
chosen temperatures over the entire frequency range measured in the 

20 experiment. A complete agreement has been found proving beyond 
reasonable doubt the validity of the present invention. 

The second set of experimental tests has been performed with pure, 
polycrystalline silicon of high purity, without the prior knowledge 
25 of the type of the majority, mobile electrical charge carriers. 
Evaporated, thick gold film or conductive silver paint were used as 
metal electrode material. 

The typical raw experimental results of the electrical response in 
30 these samples are shown in Figs. 7a and 7b, where log {Re [7 (a)]} 

S 

and log {Im[Y s (w)]} are plotted versus log of the measurement 
frequency for a number of different temperatures. 

Using the present invention it was possible to identify the 
35 electrical topology of silicon grains in polycrystalline silicon 
samples (Fig. 8a), the nature of the electrical transport in 
polycrystalline silicon, the electrical characteristics of the 
individual grains/grain boundaries (Figs. 8b, Be) and the 
concentration N of the majority mobile charge carriers (electrical 



WO 94/14078 



PCT7DK93/00410 



37 

purity) in these polycrystalline silicon samples. 

With N as free parameter an extremely good numerical fit to the raw 

experimental data (Figs. 7a ,7b) at different temperatures was 

5 obtained, using the static R,C electrical network analysis (Figs. 

8a, 8b, 8c) of the present invention. The best numerical value for the 

18 -3 

quantity N was found to be approximately 6*10 m indicating an 
extreme purity of the starting polycrystaline silicon material. 

10 Comparison between the new EIS experimental method according to the 
invention and the prior C-V (Ouasistatic. d.c. voltage bias 
dependent Capacitance measurement) and DLTS (Deep Level Transient 

spectrpscqpY) experimental methpd$ r 

15 At this stage it is appropriate to compare at least qualitatively 
the new EIS experimental method according to the present invention 
with two other experimental methods and variations thereof that are 
used today for the electrical characterization of semiconducting 
materials. It is claimed that they are capable of yielding the 

20 values of N and/or N(x) (the density of the majority mobile charge 
carriers; C-V experimental methods) and the energy position and the 
density of the deep localized energy levels within the gap of a 
measured semiconductor (the quantities N + j LQC (E^, x) and ITj LQC 
(Ej, x) in equation system (1); DLTS experimental methods). 

25 

Both of these methods attempt to measure the depletion capacitance 
C n (equation (5) and Figure (5) ) and its changes either with an 
externally applied d.c. voltage bias V? * (C-V methods) or with 
temperature T (DLTS methods). 

30 

However, they both suffer from a fundamental physical flaw in that 
the measuring frequency & at which the wanted capacitance C D of the 
system under investigation is to be measured is chosen arbitrarily. 
A quick inspection of Figure 5 shows that if this arbitrarily chosen 
35 frequency used by the two methods does not fall within a narrow 
range around fp (flat part of Cjgj ) both methods will yield incorrect 
results (the measured capacitance not related to C Q ) with deduced 
electrical material parameters that bear no relation to their true 
values. This is particularly so for the higher resistivity 
semiconductors and insulators where the usual frequency used (1 MHz) 
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would yield a measurement of the bulk capacitance- Cg rather than 
that of the depletion capacitance C Q (equation (5) and Figure 5). 

It is also evident that when the measuring frequency is chosen 
correctly (close to the flat part of the - Figure 5), both 
methods are represented by a single point on the measured electrical 
impedance curve Z (a>) (Im(Y(a>)/u)) of the present invention. In such 
case they therefore represent a measurement of the level (Figure 
5) as a function of the d.c. voltage bias and/or as a function of 
the temperature T. 

The knowledge and the proper analysis of the full Z (w) curve 
(present invention) is however essential in order to decide about 
the validity of the C-V and DLTS results. 

It should also be pointed out that the EIS experimental method 
according to the present invention measures N and/or M(x) directly 
through the measurement of the I (u) over the entire relevant 
frequency range. One of the major results of the present invention 
1s namely the discovery of one to one correspondance between the 
measuring frequency u and the spatial co-ordinate x. The measurement 
of Z (w) at frequency ^ Is a measure of the density of the mobile 
electrical charges N(x) at point x^. This is a result of the solu- 
tion of the Poisson equation and it 1s valid under small signal 
approximation - the usual condition of the EIS experiment. 

In conclusion it should be stressed again that the C-V and DLTS type 
of measurements are embodied automatically (as a single point in the 
Z (<a) curves) in the EIS experimental method according to the 
30 present invention. No dedicated C-V or DLTS apparatus is necessary 
and by measuring T (a) in the whole frequency range there is a 
quarantee that no errors arise from incorrectly chosen measuring 
frequency. 

35 References: 
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Patent Claims. 

1. A method for determining characteristic electrical properties of 
semi -conducting materials wherein the time or frequency dependent 
5 electrical impedance or admittance of the material is measured, 
characterized by the steps of establishing a measurable 
difference in chemical potential between the inner and surface of 
the material or between the inner of the material and the metal 
electrode - material surface interface layer, and of determining 

10 from the measured electrical impedance or admittance one or more of 
the following electrical parameters of the material: 
the difference A/i ch in chemical potential, 
the dielectric constant £, 
the density N of majority mobile charge carriers, 

15 the density N - of minority mobile charge carriers, 

the electrical mobility /* of majority mobile charge carriers, 
the electrical mobility ^ n of minority mobile charge carriers, 
the emission and capture rates E.R and C.R, respectively, for mobile 
positive and negative charge carriers characterizing the effect of 

20 surface and bulk localized states within the band gap (Eg), when 
they are present, on the electrical transport and the measured 
electrical impedance. 

2. The method according to claim 1, characterized in 
25 that the determination of the electrical parameters of the material 
is carried out by solving a system of equations for the total charge 
density p(x,t) consisting of the mobile negative and positive charge 
densities and localized negative and positive charge densities 1n 
the material supplemented with initial and boundary conditions, 
30 where the space-time development of p(x,t) determines the electric 
current running in the external circuit in response to the applied 
electric voltage thereby defining the electrical complex impedance 
Z $ (w) of the material. 

35 3. The method according to claims 1 and 2, characterized 
in that the electrical parameters of the metal electrode-material 
surface interface layer such as the emission and capture cross 
sections a e and a c of the Interface layer for the mobile positive 
and negative charge carriers are determined from the measured 
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electrical Impedance by means of the initial and boundary conditions 
for p(x,t). 

4. The method according to claims 1 and 2, in which the interface 
layer is a Si -SiO,, -metal electrode interface layer, charac- 
terized in that the electrical elements and of the 
interface layer are determined by means of the expressions: 



e P l- ( cy*T) 2 m ± m sinQt-cyJcT) . ^ 
Cl -\e\-A*T 2 A r (K-cyJrtT) 



15 C - e • ^ 



r 



where 



\ - ~ ' . « - 4' *• and ^ = « • d sio , 



20 Cl = T 



The quantities appearing in the analytical expressions 
25 have the following meaning : 



c sio? - dielectric constant of the silicon dioxide 

d si0 ^ - thickness of the native silicon dioxide layer 

30 A r - total area of the sample 

A m Richardson constant 

AE X - quantum tunnelling barrier height 

m - electron rest mass 

h - Planck's constant 

k - Boltzraann's constant 

35 T - absolute temperature 

AE 2 - activation energy 
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5. The method according to claim 2, characterized in 
that the electrical parameters of the material are determined from 
the following expressions for capacitance and resistance of the 
inner of the material and the depletion region thereof, 
respectively: 



10 



Rb |e| V (t) -n 



L - 2'WM 



c * ~ e °"** " L - 2-w D cn 



15 



Kp - 



2-k-T 



C D = e o*«r 



\e\y/ir e 'N v v.iTfVblTI 



k-T 
2^11^(7) 



arctan 



W D (T) 



20 



25 



30 

which are identifiable as levels on the measured electrical 
impedance curve Z («). (Fig. 5 and equation (3)). 

35 

(The meaning of the individual quantities in these formulas are 
explained elsewhere in the text). 
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6. The method according to claim 2, characterized in 
that the electrical parameters of the material are determined from 
the capacitance and resistance of the Interface layer and the Inner 
of the material and the electrical impedance of the depletion region 
5 by the following expressions: 



10 



*-H+S 



c^e^A'-T 2 A z (ncy/c-T) 



y -M*S t SiO, ^ 



15 



^ 1 + 0) 2 Ti 



arc tan 



1 ♦ OJ 2 T* 



c ( 1 + i * ox) j 



20 



B |e|Ti(r>-AT A r 



L - 2-W D (T) 



25 where the values of and R M+S are identifiable as levels on the 
measured impedance curve Z(w) if such levels are present. 

(The meaning of the individual terms in these expressions are 
explained elsewhere in the text). 

30 

7. The method according to claim 1, characterized in 
that for the "research" experimental set-up the material is 
deposited semipermanently on the metal electrodes or vice versa. 

35 8. The method according to claim l f characterized in 
that for "quality control" set-up the metal electrodes are brought 
into contact with the surface of the material via an electrically 
conducting medium which ensures a well-defined electrical contact 
area on the specimen, such as an elastic or flowable material which 
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may be brought into and out of engagement with one or both surfaces 
of the material . 

9. The method according to claim 8, characterized in 
5 that the measurement is repeated in points distributed throughout 

the area of the specimen for mapping of the electrical parameters of 
the material. 

10. An apparatus for carrying out the method according to claim 1, 
10 comprising an external electrical circuit for measuring the time or 

frequency dependent electrical impedance or admittance of the 
material and means for electrically contacting the surface or 
surfaces of the material, characterized by means for 
providing a measurable difference in chemical potential between the 
15 inner of the material and its surface or between the inner of the 
material and the metal electrode - material surface interface layer, 
and means for determining from the measured electrical impedance or 
admittance one or more of the following electrical parameters of the 
material : 

ch 

20 the difference A/r in chemical potential, 
the dielectric constant £, 

the density N of majority mobile charge carriers, 

the density N m j n of minority mobile charge carriers, 

the electrical mobility n of majority mobile charge carriers, 

25 the electrical mobility of minority mobile charge carriers, 

the emission and capture rates E.R and C.R, respectively, of mobile 
positive and negative charge carriers characterizing the effect of 
surface and bulk localized states within the band gap, when they are 
present, on the electrical transport and the measured electrical 

30 impedance. 

11. The apparatus according to claim 10, characterized 

in that it is arranged for measuring the electrical 

Impedance/admittance 1n a frequency range from 1 mHz to 10 GHz 

35 (1000 s - 0.10 ns) and with a real component of the Impedance in the 

14 

range from 0,001 ohm to 10 ohm and an Imaginary component of the 
admittance divided by the frequency o> in the range from 10" 14 C to 

io- 2 . 
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12. The apparatus according to claim 10, characterized 
by comprising a specimen support having a plurality of coaxial leads 
for connection to the specimen of the external electrical 
measurement circuit and means for controlling the temperature of the 

5 specimen. 

13. The apparatus according to claim 10, characterized 
by comprising electrically conducting contact means providing for a 
well defined electrical contact area on the specimen and which are 

10 arranged for being brought into and out of engagement with one or 
both of the surfaces of the material, and control means for 
relatively adjusting the contact means and the specimen throughout 
the total area thereof. 

15 



20 



25 



30 



35 
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